Direct Observation of the Hyperfine Transition of the Ground State Positronium 



T. Yamazaki, 1 A. Miyazaki, 1 T. Suchara, 1 T. Namba, 1 S. Asai, 1 T. 

Kobayashi, 1 H. Saito, 2 I. Ogawa, 3 T. Idchara, 3 and S. Sabchcvski 4 
1 Department of Physics, Graduate School of Science, 
and International Center for Elementary Particle Physics, 
University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

2 Institute of Physics, Graduate School of Arts and Sciences, 
University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan 

^Research Center for Development of Far-Infrared Region, 
University of Fukui, 3-9-1 Bunkyo, Fukui-shi, Fukui 910-8507, Japan and 
4 Bulgarian Academy of Sciences, 72 Tzarigradsko Shose Blvd., 1784 Sofia, Bulgaria 

(Dated: March 2, 2013) 

We report the first direct measurement of the hyperfine transition of the ground state positronium. 
The hyperfine structure between ortho-positronium and para-positronium is about 203 GHz. We 
develop a new optical system to accumulate about 10 kW power using a gyrotron, a mode converter, 
and a Fabry-Perot cavity. The hyperfine transition has been observed with a significance of 5.4 
standard deviations. The transition probability is measured to be A = 3.llj 2 x 10 -8 s _1 for the 
first time, which is in good agreement with the theoretical value of 3.37 x 10~ 8 s _1 . 



Positronium (Ps) [l[ , a bound state of an electron and 
a positron, is a purely leptonic system and is a good tar- 
get to study quantum electrodynamics (QED) in bound 
state. The triplet (l 3 Si) state of Ps, ortho-positronium 
(o-Ps), decays into three gamma rays with a lifetime of 
t = 142 ns [1,11. On the other hand, the singlet (l 1 ^) 
state of Ps, para-positronium (p-Ps), decays into two 
gamma rays in r p = 125 ps Q. The energy level of the 
ground state o-Ps is higher than that of the ground state 
p-Ps due to the spin-spin interaction between the electron 
and the positron. This difference is called the hyper- 
fine structure of the ground state positronium (Ps-HFS), 
which is about 203 GHz. Although precise measurements 
of Ps-HFS have been performed in 1970s and 1980s [1,0], 
all of them are indirect measurements using Zeeman split- 
ting of about 3 GHz caused by a static magnetic field of 
about 1 T. There is a discrepancy of 3.9 standard devia- 
tions (15 ppm) between the measured and the theoretical 
value [Tj . The largest systematic uncertainty common to 
all previous measurements is the non-uniformity of the 
static magnetic field. It is important to directly measure 
Ps-HFS, in order to avoid the systematic uncertainty of 
the static magnetic field. Here we present a direct obser- 
vation of the hyperfine transition between Ps-HFS, which 
is the first great step toward a direct measurement of Ps- 
HFS. The hyperfine transition of the ground state Ps, 
which is Ml transition, has not yet been observed di- 
rectly, since the transition probability (Einstein's A co- 



efficient is A = 3.37 x 10" 



is 10 times smaller 



than the decay rate of o-Ps (7.0401(6) x 10 6 s" 



0,1). In 



order to cause sufficient amount of stimulated emission 
from o-Ps to p-Ps, we develop a new optical system which 
consists of a gyrotron as a sub-THz radiation source, 
a mode converter to convert the gyrotron output to a 
Gaussian beam, and a Fabry-Perot cavity to accumulate 
high power sub-THz radiation. The gyrotron is a novel 
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FIG. 1: Schematic diagrams of our experimental setup. Top 
view of the gas chamber is shown in the box. Ml and M2 
are parabolic mirrors made of aluminum. We use a gold mesh 
plane mirror with a transmittance of about 3 % as a beam 
splitter (BS). Three pyroeiectric detectors (PY) are used to 
monitor the incident, the reflected and the transmitted power. 



high power radiation source for sub-THz to THz region, 
which enables us to perform a direct measurement of the 
hyperfine transition. High power 203 GHz radiation in 
the Fabry-Perot cavity causes the hyperfine transition 
from the ground state o-Ps to p-Ps, and p-Ps promptly 
decays into two back-to-back 511 keV gamma rays. Con- 
sequently, the transition signal (o-Ps — > p-Ps — > 2j) has 
distinctive features that it has a lifetime of o-Ps and de- 
cays into two back-to-back 511 keV gamma rays as p-Ps. 

Figure [TJ shows a schematic view of our experimental 
setup. We use Gyrotron FU CW V which produces 
202.89 GHz (140.06 GHz) radiation in TE 03 (TE 02 ) mode 
in 15 ms pulses at 20 Hz. The power is monitored with a 
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pyroelectric detector, which is fed back to voltage of the 
heater of the electron gun. As a result, it can operate 
stably with about 300 W power within 10 % fluctuation. 

In order to enhance the output power of the gyrotron, 
the radiation is accumulated in a Fabry-Pcrot cavity. The 
gyrotron output (TEo„ mode) is converted to a Gaussian 
beam so as to obtain good coupling with the Fabry-Pcrot 
cavity. Main components of a mode converter are a step- 
cut waveguide and a large parabolic mirror made of alu- 
minum (Vlasov antenna). They convert TEo„ mode to 
a bi- Gaussian beam geometrically if the axis of the step- 
cut wave guid e and the focal point of the parabola are 
matched U0JI. Two mirrors (Ml and M2) are used to 
convert the bi-Gaussian beam into a Gaussian beam. In 
order to improve the beam quality, we insert an aper- 
ture (diam. = 50 mm) as a spatial filter to block out side 
lobes of the beam. Spatial distribution of the beam is 
measured by exposing a PVC sheet to the beam and tak- 
ing its picture by an infrared camera. Power conversion 
efficiency, which is estimated from the spatial distribu- 
tion, is 28 ± 2 % due to a limitation of a purity of wave 
mode in the gyrotron output. 

The Fabry-Perot cavity is made with a gold mesh plane 
mirror (diam. = 50 mm) and a copper concave mirror 
(diam. = 50 mm, curvature = 300 mm). The incident 
Gaussian beam resonates within the Fabry-Perot cav- 
ity when the cavity length (136 mm) is equal to a half- 
integer multiple of the wavelength of the radiation (about 
1.5 mm). The cavity length is controlled by moving 
the copper concave mirror mounted on an X-axis stage 
(NANO CONTROL TS102-G). The gold mesh plane mir- 
ror is a key component of the Fabry-Pcrot cavity, and is 
made on a Si02 plate using photolithography and liftoff 
technique. The line width and separation are 200 /im 
and 160 fim, respectively. The mesh parameters are de- 
signed to obtain high reflectivity (99.38 %) and reason- 
able transmittance (0.39 %), which are simulated with 
CST Microwave Studio [Hj]. As a result, the finesse of 
the Fabry-Perot cavity attains T = 623 ± 29, which is 
estimated from the width of the resonance peak while 
changing cavity length. The power accumulated in the 
Fabry-Perot cavity reaches about 10 kW. 

The power accumulated in the Fabry-Perot cavity is 
estimated with the power transmitted through a hole 
(diam. = 0.6 mm) at the center of the copper con- 
cave mirror. The transmitted power is monitored with 
a pyroelectric detector. A ratio between accumulated 
power and transmitted power is obtained from indepen- 
dent measurements using the Gaussian beam as follows. 
First, the beam undergoes total absorption in water. Its 
total power is estimated from a temperature increase of 
the water. Next, the copper concave mirror is exposed 
to the beam, and power transmitted through the hole 
is measured with the pyroelectric detector. From these 
measurements, the ratio of the Gaussian beam power to 
transmitted power is obtained. However, the spatial dis- 



tribution of the Gaussian beam is different from that of 
the beam inside the Fabry-Perot cavity. Correcting the 
difference of the beam shapes and considering that only 
the beam going to the copper mirror direction can be 
transmitted through the hole, the ratio of accumulated 
power to transmitted power is obtained. The Gaussian 
beam shape is measured with a PVC sheet, and the spa- 
tial distribution in the Fabry-Perot cavity is calculated 
from the cavity length, the curvature of the copper mir- 
ror, and the radiation wavelength. The uncertainty of 
the ratio of accumulated power to transmitted power is 
-30 % because of the fluctuation of the beam shape be- 
tween the exposure on the copper mirror and that on 
the PVC sheet. The uncertainty of power does not affect 
the direct measurement of the hyperfinc transition, but 
contributes to the accuracy of the transition probability, 
which is also measured in our experiment. 

Positronium formation assembly shown in Fig. [1] is as 
follows: A 780 kBq 22 Na positron source is placed above 
a thin plastic scintillator (NE-102, thickness = 0.1 mm). 
Emitted positrons pass through the scintillator and pro- 
duce light pulses that are directed to two 1.5-inch fine- 
mesh photomultipliers (Hamamatsu R5924-70) by the 
light guide. Positrons form Ps when stopped in the mixed 
gas (1.9 atm N 2 and 0.1 atm i-C 4 Hi ) [12]. About 5 % 
of positrons are tagged by the plastic scintillator and 
stop in the gas, and then, about 1/4 of them form Ps. 
Therefore, the Ps formation rate is about 10 4 s _1 . Ps 
has kinematic energy of about 1 eV just after its forma- 
tion. It becomes thermalized after O(10 ns) with clastic 
collisions with gas molecules and the kinetic energy be- 
comes about 1/30 eV. Since we use delayed coincidence 
as shown in Fig. [2l the width of the Doppler broad- 
ening due to motion of thermalized Ps is only about 
A/d = 0.08 GHz, which is much smaller than the natural 
linewidth A/ n - 1/2ttt p = 1.27 GHz. 

Gamma rays emitted from Ps decay are observed in 
four LaBr 3 (Ce) crystals (Saint-Gobain Crystals, diam. = 
1.5 inch and length = 2.0 inch). The four detectors are 
placed as shown in Fig. [1] to make four back-to-back 
pairs. The scintillation pulses of the LaBr 3 (Ce) crys- 
tals are detected with 1.5-inch fine-mesh photomultipli- 
ers (Hamamatsu R5924-70). The energy resolution of the 
LaBr 3 (Ce) detectors is 4 % (FWHM) at 511 keV. The 
primary decay time is 16 ns. These are advantages for 
tagging monochromatic 511 keV gamma rays and avoid- 
ing pileup of gamma ray signals. 

In order to select Ps decay events, data acquisition 
logic is set up as follows: When at least one back-to-back 
signal from the LaBr 3 (Ce) scintillator pairs is coincident 
within 40 ns, and then when this coincidence is within 
— 100 ns to 1100 ns of the timing of the plastic scintillator, 
data acquisition is triggered. A charge ADC (PHILLIPS 
7167) and another charge ADC (REPIC RPC-022) are 
used to measure the energy information of the plastic 
scintillator with short and long gate, respectively. The 
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energy difference between short and long gates is used to 
suppress accidental background as mentioned later. The 
outputs of the LaBr 3 (Ce) detectors are recorded with a 
charge ADC (CAEN C1205). The time information be- 
tween the plastic and LaBr3(Cc) scintillators is recorded 
using a direct clock (2 GHz) count type TDC (KEK 
GNC-060) 0. 

Four runs have been performed. In three runs (Run I, 
III, and IV), 202.89 GHz radiation (TE 03 mode) is used, 
and different powers are accumulated in the Fabry-Perot 
cavity (11.0 kW, 0.0 kW, and 5.6 kW). In another run 
(Run II), off-resonance frequency of 140.06 GHz in TE02 
mode is used to check systematic effects due to the ab- 
sorption of the radiation in the mixed gas. Total period 
of data acquisition is about two weeks. During the data 
acquisition, energy and time calibrations arc performed 
every 30 minutes. Trigger rates are about 1 kHz. The 
7-ray peak at 511 keV and the zero energy peak are used 
to calibrate the LaBr3(Ce) detectors. The room temper- 
ature is maintained within 26 ±1 °C in order to maintain 
good stability during the data acquisition. 

Figure [2] shows the time difference between the plas- 
tic scintillator signal and the coincidence signal of the 
LaBr3(Ce) detectors. A sharp peak from prompt anni- 
hilation is followed by the exponential curve of transi- 
tion signals and o-Ps decay signals, and then the con- 
stant spectrum due to accidental overlaps of a triggered 
positron and uncorrelated gamma rays. A good timing 
resolution (a = 0.8 ns) is obtained. After selecting a time 
window from 50 ns to 350 ns to enhance the transition 
signals and o-Ps decay events, accidental events remain 
as the dominant source of back-to-back 511 keV gamma 
rays. In the case of accidental events, there is another 
plastic scintillator hit at the timing of 7-ray hit. The 
energy deposit on the plastic scintillator measured with 
long gate becomes larger than that measured with short 
gate. To reject accidental events, the energy difference 
between long gate and short gate is limited from —2.5 pe 
(photoelectron) to 1.7 pe. This cut is applied on both 
PMTs of the plastic scintillators. 

Finally, we count the number of events in which back- 
to-back 511 keV gamma rays are observed. Figure [3] 
shows the energy spectra measured with the LaBr3(Ce) 
scintillator in the highest power on-resonance run (Run 
I, II.O+33 kW). The delayed coincidence and the ac- 
cidental rejection arc applied. In addition, a 511 keV 
7-ray hit on the LaBr3(Ce) scintillator at the opposite 
side of the back-to-back pair is required, where the en- 
ergy window is set from 494 keV to 536 keV. Remain- 
ing accidental background is estimated from the events 
in another time window set from 850 ns to 900 ns, and 
is subtracted. Circles and triangles show "beam ON" 
and "beam OFF" spectra, respectively. The data taken 
during "beam OFF" period in the pulse beam are used 
to estimate background. The "beam OFF" spectrum 
consists of pick-off annihilation (o-Ps + e~ — > 2j + e~) 
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FIG. 2: Time difference between the plastic scintillator and 
the coincidence signal of the LaBrs(Ce) scintillators. Solid 
line and hatched histogram show the time spectrum before 
and after accidental rejection cut, respectively. The time win- 
dow for delayed coincidence is shown as a dashed line. 
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FIG. 3: Energy spectra of the LaBra(Ce) scintillator in the 
highest power on-resonance run (Run I, 11.0i|'g kW) after 
all event selections are applied. Circles and triangles show 
"beam ON" and "beam OFF" spectra, respectively. 



and 37 decay (o-Ps — > 37) of o-Ps. Transition signals 
(o-Ps — > p-Ps — > 27) increase when o-Ps are exposed to 
high power sub-THz radiation during "beam ON" period. 
The signal rate in the energy window from 494 keV to 
536 keV is i?oN — Roff = 15.1 ± 2.7(stat.) mHz, where 
-Ron (#off) is the "beam ON" ("beam OFF") event rate 
after all event selections are applied. 

Systematic errors are summarized in Table HI The 
largest contribution is the uncertainty in Ps formation 
probability. Ps formation probabilities of the "beam ON" 
and the "beam OFF" data are different because of ab- 
sorption of the sub-THz radiation in the mixed gas, which 
is enhanced when the beam resonates with the Fabry- 
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TABLE I: Summary of the systematic errors. The values are 
ratios to the background. 



Source 


Run I 


Run II 


Run III 


Run IV 


Ps formation prob. 


-0.27 


% 


-0.39 


% 


+0.20 


% 


-0.13 


% 


Accidental rejection 


+0.17 


% 


+0.05 


% 


+0.13 


% 


±0.23 


% 


Energy resolution 


-0.08 


% 


+0.06 


% 


-0.11 


% 


-0.02 


% 


BG normalization 


±0.03 


% 


±0.04 


% 


±0.04 


% 


±0.04 


% 


Total 


+0.17 

-0.29 


% 


+0.08 
-0.39 


% 


+0.24 
-0.12 


% 


+0.24 
-0.14 


% 



Perot cavity. The difference is estimated by counting the 
number of events in the time window before the energy 
cut is applied, since Ps formation probability is indepen- 
dent of the 7-ray energy cut. The difference in Ps forma- 
tion probability is the largest in off-resonance run (Run 
II). Another dominant systematic error is uncertainty in 
the efficiency of the accidental rejection cut. Inefficiency 
of the accidental rejection depends on the rates of the 
plastic scintillator signals which go over the discriminator 
threshold (~ 1 pe). This systematic effect is estimated 
from the difference of the efficiency of the accidental re- 
jection between "beam ON" and "beam OFF" , which 
is independent of the 7-ray energy cut. In addition, if 
the energy resolution and energy scale of the LaBr 3 (Ce) 
scintillator are different between "beam ON" and "beam 
OFF", fake signals appear because of the back-to-back 
511 keV energy selection. This effect is estimated using 
Monte Carlo simulation with Geant4 where the en- 
ergy resolution and energy scale taken from data are used 
as input. The last dominant source is the uncertainty of 
background normalization. The background is estimated 
from "beam OFF" events. Its normalization is performed 
using the number of events in the prompt time window 
set from —3 ns to 1.5 ns, where the usual e + annihilation 
is dominant (77 %). Statistical accuracy determines the 
normalization uncertainty. 

The systematic errors discussed above are indepen- 
dent, and the total systematic error can be calculated as 
their quadrature sum. Final result with the systematic 
errors is 



Ron — Roff 



15.1±2.7(stat.)±g;|(sys.) mHz. (1) 



This is the first direct observation of the hypcrfme tran- 
sition of the ground state positronium with a significance 
of 5.4 standard deviations. In addition, the fraction of 
the transition signals is proportional to the power accu- 
mulated in the Fabry-Perot cavity (Fig. 0J, and the off- 
resonance data (Run II) gives a null result as expected, 
despite the relatively large difference in Ps formation 
probability as seen in Table |TJ 

The transition probability (or Einstein's A coefficient) 
between the ground state Ps-HFS is also measured for the 
first time. It can be estimated from the observed transi- 
tion rate, the power accumulated in the Fabry-Perot cav- 
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FIG. 4: Power dependence of the amount of the transition sig- 
nals. The vertical axis shows signal to background ratio. The 
horizontal axis shows the power accumulated in the Fabry- 
Perot cavity. The error bars contain statistical uncertainty as 
well as systematic uncertainties summarized in Table |U The 
dashed line shows the result of a linear fit. 



ity, and 27/37 detection efficiency estimated from Monte 
Carlo Simulation with Geant4. The estimated result is 



A = 3.1 



+1.6 



x 10" 



(2) 



which is consistent with the theoretical value of 3.37 x 
10 -8 s _1 [1]. The largest uncertainty is the estimation 
of the absolute power accumulated in the Fabry-Perot 
cavity. 

Our next target is to directly measure Ps-HFS for the 
first time. Output frequency of gyrotron can be changed 
with cavities of different sizes. In Ps-HFS measurement, 
relative accuracy of the power estimation at different 
frequency points is necessary. In addition, in order to 
perform precise measurement of Ps-HFS, we need more 
statistics. A possible way to increase statistics is to use 
a slow positron beam and make positroniums in vacuum 
using a thin metal foil [l4[ . It also eliminates systematic 
uncertainty and beam power loss due to absorption of 
the sub-THz radiation. 

In summary, the hypcrfme transition of the ground 
state positronium has been observed directly for the first 
time with a significance of 5.4 standard deviations. We 
develop a new optical system to accumulate about 10 kW 
power using a gyrotron, a mode converter, and a Fabry- 
Perot cavity, in order to cause observable amount of stim- 
ulated emission from o-Ps to p-Ps. The transition prob- 
ability (or Einstein's A coefficient) is also measured to 
be A = 3.lj^ 2 x 10~ 8 s _1 for the first time, which is in 
good agreement with the theoretical value. 

This research was funded in part by the Japan Society 
for the Promotion of Science. 
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